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Abstract 

In metallic nanoparticles, shifts in the ionization energy of surface atoms with 
respect to bulk atoms can lead to surface bands. Within a simple Tight Binding model 
we find that the projection of the electronic density of states on these sites presents two 
overlapping structures. One of them is characterized by the level spacing coming from 
bulk states and the other arises from the surface states. In very small particles, this 
effect contributes to an over-broadening of the NMR absorption spectra, determined 
by the Knight shift distribution of magnetic nuclei. We compare our calculated Knight 
shifts with experiments on aluminum nanoparticles, and show that the deviation of the 
scaling law as a function of temperature and particle size can be explained in terms of 
surface states. 
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1 Introduction 



Several NMR experiments on metallic nanoparticles have shown PjJ |2l EI quantum size ef- 
fects: the absorption spectra becomes quite broad and asymmetric when either temperature 
or particle size diminishes. Studies of relaxation time |H| T\ indicate that the spectra is 
inhomogeneous, implying phenomena of a local nature. This inhomogeneity has been inter- 
preted jl] considering a strictly local description for the electron polarization in the presence 
of a magnetic field H. Given the Bohr magneton, /ib, and N(ef,t), the local density of s 
states at the Fermi level e-p, the local Pauli susceptibility is written as: 



This equation expresses that, due to the finite size effects and interference phenomena proper 
of the mesoscopic systems, the electronic spin polarization is inhomogeneous. The hyperfine 
coupling between nuclear spins and electronic spins produces different shifts in the resonance 
frequency for each nucleus in the nanoparticle: 



Therefore, an inhomogeneously broadened NMR absorption line is the key to indirectly 
measure the fluctuations of the local density of states (LDOS). 

In previous works [S] we have calculated the NMR line shape in metallic nanoparticles 
within a Tight Binding model with M orbitals accounting for the fluctuations of the LDOS. 
A remarkable outcome, consistent with many experimental results [DE], is a universal scaling 
behavior of the line position and shape with respect to the variation of the thermal energy 
kftT and the mean spacing level A = 1/ [M x Nq^sf)]. This results in a universal scaling 
law for the Knight shifts and the Pauli susceptibility 6J whose relevant parameter is 






(2) 



a = k B T/A. 



(3) 
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However, in very small crystalline particles, the experimental line shape presents an anoma- 
lous over- broadening at intermediate temperatures [7j . In addition, a.c. conductivity mea- 
surements show that tunneling among particles plays a relevant role in the electronic prop- 
erties |8j. This suggests the need to explore for additional sources of fluctuations on the 
statistics of energy levels, particularly the effects of surface states. It is known that self 
consistent calculations of metal surfaces [3 E] give different ionization energies for surface 
and bulk sites. This fact is expected both in metallic particles with clean surfaces and as 
well as those with chemisorbed atoms. In our previous calculation we made no attempt to 
consider surface states. However the natural connectivity of the surface orbitals caused slight 
departures from the universal scaling law [5]. 

In the present work we use a Tight Binding Hamiltonian with an ad hoc shift in the surface 
site energies [TUHH]. This produces surface states, and as a consequence, large fluctuations on 
the LDOS distribution. According to the Eqs. Q and El this leads to a modified distribution 
of Knight shifts detected as a wider NMR absorption lines. We compare our calculations of 
the Knight shift as a function of the scaling parameter with experimental values on aluminum 
nanoparticles and show that the observed deviation of the scaling law is a manifestation of 
surface effects. 

2 Model Hamiltonian 

A metal particle with M atoms can be modeled with the Hamiltonian 

M M M 

n = Y, E *4ci + Yl + v^tcj), (4) 

i=l j>i i=l 

where Ei is the energy of an s state centered at site i of a cubic lattice. Vij = V is the kinetic 
energy involved in hopping between nearest neighbors sites i and j. To represent shape and 
crystal inhomogeneities, the sites energies are taken in the range — W/2 and W/2 (Anderson's 
disorder). We define surfaces sites as the sites which being at the surface of the cube, have an 
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additional energy shift U. The rest of the sites are bulk sites, even though they may belong to 
the faces of the cube with un-shifted site energies. This allows to model different surface to 
bulk ratios maintaining M constant. Tunneling among particles produces an inhomogeneous 
broadening [IS] of the atomic energy levels (Ei — » E^— iTj), taking its higher values T s at 
particular sites (contacts) at the surface. Assuming that Tj <C V, and neglecting localization 
effects we can assume Ij = r] « T s /M for all sites in the nanoparticle. As shown in the Ref. 
4 finite temperature effects are included by an additional broadening on the energy levels, 
V — Vo + k&T. Properties of the single particle excitation spectrum are contained in the 
retarded (advanced) Green's function 



G 



R(A), 



E 



s + [(-)i7J - E k 



(5) 



where ipk( T ) — '^2i a k( r i) ( Pi( r ) an d E^ are the exact eigenf unctions (molecular orbitals) and 
eigenenergies for the isolated particle, respectively, rj is a natural broadening of the electronic 
states and the (— ) sign corresponds to the retarded Green's function. The local density of 
states per atom (LDOS) at the i-ih site is evaluated as: 



N(e,Ti) 



2^[GUe)-GUe) 



(6) 



from which the relevant contribution to the density of states per unit volume at the i- 
th nucleus, iV(e, rj),is obtained. The evaluation of the local Green's function in Eq. |S] is 
obtained via the Matrix Continued Fraction method JH]- Its basic idea is to exploit the 
short range interactions in the Hamiltonian (@J) by indexing states in a way that subspaces 
representing layers interact through nearest neighbor subspaces. In matrix form: 
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where O's are null matrices, E's in the diagonal represent intra-layer interactions while the 
only non-zero off-diagonal matrices "V n ,n±i connect nearest neighbor layers. Detailed struc- 
ture of the sub-matrices depends on the lattice, for the cubic structure V Hjn ±i = VI, with 
1 the identity matrix. The local retarded Green's functions connecting sites % and j within 
the n-th layers are arranged in a matrix 

where the matrix self energies and are calculated in terms of Matrix Continued 
Fractions (MCF) defined through the recurrence relations: 

± = V n , n ±l - R ±V n ±i )Tl , (9) 

(£ + lf]j 1 — H(n±l,n±l - ^n±l 

which are calculated with the boundary conditions: = = 0, where L denotes the 
number of layers. 

Model parameters. Parameters that model an s-band are ^3]: V = 0.9eV, which is 
consistent with a bandwidth of B = 12V = 0.8Ry = lleV; a shift U = 2V, consistent with 
a shift of 0.12Ry estimated for metal clusters and a Fermi energy of ep = 3V + B/2, which 
measured from the band bottom gives a ratio e-p/B = 0.75, close to typical values of metals 
[14J. In this work we use rj = 0.051^, which assumes that the main contribution to the level 
broadening comes from tunneling i.e. r] > k^T which is consistent with the conductivity 
measurements [S]. 

The statistical distribution of the LDOS is calculated taking an ensemble of ten Fermi 
energies in a range Sep = 0.5V around £p. For each energy ten disorder configurations are 
considered. Each LDOS N(ep, r) is normalized to the bulk value N (ep), which is evaluated 
at the central site of a particle with 15 x 15 x 15 orbitals. The normalized LDOS is defined 
as x = N(ep,r)/N (ep), and occurs with a probability I(x). According with Eqs. [T]and|21 
I(x) is also the absorption at normalized frequency shifts: x = {oj — u )/(u K — u) ). Here 
uj k is the bulk metal frequency and uj the frequency of metal nuclei in dielectric materials 



(i.e. salt). Due to this correspondence, from now on, I(x) will refer indistinctly to either the 
distribution of the normalized LDOS at the Fermi energy or the Knight shift NMR spectrum. 
Quantum size effects manifest as shifts of the NMR line maximum: as the particle size or 
temperature decreases, the shift goes from x = 1 (metal bulk) to x = (salt) while the line 
broadens and turns asymmetrical jHEl- It is important to emphasize that the simple model 
used in this work only attempts to provide physical insight into the quantum size effects as 
manifested in the LDOS distribution, as well as a qualitative description of surface effects. 

3 Results 

To model small surface to volume ratios proper of a big particle, the surface sites (with mean 
site energy U = 2V) were chosen on a single face of a 7x 7x 7 cube, i.e. M sur f = 7x 7 = 49 are 
surface sites and Mbuik = 6x7x7 = 294 are bulk sites. In this case the ratio of level spacing 
is A sur f./Abuik = 4. Results for this particle can be extrapolated to an actual particle with 
40 x 40 x 40 orbitals having roughly the same surface to bulk ratio. Figure 1 shows the Knight 
shift spectrum obtained from the LDOS occurrence distribution at surface sites (thick line). 
It shows that the shift in the site energies produces an over-broadening with respect to the 
Knight shift distribution from the same sites for the homogeneous configuration with U = 
(thin line). The over-broadening suggests the existence of multiple fluctuation scales in the 
NMR spectrum. To clarify this, Figure 2(a) shows the projection of the density of states 
(DOS) at surface sites, N sur{ (e). It reveals two structures, one corresponding to surface sites 
with a typical spacing A surf . and another one with typical spacing A bulk corresponding to 
bulk sites. Notice that the structure of bulk states (Fig. 2.b), characterized by a mean level 
spacing of A bulk , enters into the surface DOS as a substructure. These bulk states have 
a small weight on the surface band because U provides a barrier that prevents the mixing 
between states with the same kinetic energy (E^,,) parallel to the surface. 

The isolated surface and bulk bands iV s urf (e) and A^ ulk (e), provide a rationale to un- 
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Figure 1: Knight shift spectrum of surface sites when U = IV (thick line), compared with 
the case with U = (thin line). M = 7 x 7 x 7, M surf = 49 (corresponding to one side of 
the cube) and M bu i k = 294. W = IV and r] = 0.05V 

derstand the structure of the surface DOS, and consequently the multiple fluctuation scales 
in the NMR spectrum. Typically A?" surf .(e) = aiV s urf (e) + bN® ulk (e). Within a perturbative 
calculation, a < 1 — (V 2 /U) 2 = 3/4 and a + b = I. The main peaks in A^ surf .(5) (Fig. 2. a) 
arise from those in N® uii (e). They produce a wide range of high LDOS values with low prob- 
ability (tail in the J surf .(a;), Fig. 1). The valleys of the global structure will produce lower 
density values with high probability (peak in the / S urf.(^))- In the same way, the denser spec- 
trum of the bulk band N® u]k (e), which contributes with small weight, provides a lower scale 
fluctuations to N Bxa f m (e). Again, the peaks and valleys of iV° ulk (e) produce an asymmetric 
distribution Jbuik(^) with a maximum. Therefore, the superposition of these two LDOS will 
cause the over-broadening in the distribution of Knight shifts for the surface sites roughly 
represented by: 




(10) 
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Figure 2: DOS projection on a) surface sites with energies U = 2V and b) bulk sites with 
energies U = 0. M = 7 x 7 x 7, M surf . = 49 and M bulk = 294. W = IV and rj = 0.05V. 



Even without attempting a fitting, this shows that the predictions of numerical results are 
consistent with the spectral analysis. 

In order to evaluate the effect of surface states on bulk NMR signal of big particles, 
Figure 3 shows the Knight shift spectrum (thick line) for the bulk sites together with a case 
with unperturbed surfaces (U = 0) represented by the thin line, ibuik(^) does not present 
much over-broadening with respect to the homogeneous configuration. The reason for this 
can be understood by analyzing the DOS projection on bulk sites. The relevant scale in 
^Vbuik(^) (Fig. 2.b) is Abuik- Since surface states enter into the bulk with small weight, and 
since A sur f. > Abuik (A sur f. = 4Ab u ik), the relevant scale is dictated by spacing of bulk levels. 
Therefore, the surface band does not affect the statistics of the LDOS of bulk sites. 

The previous analysis suggests that if the surface to volume ratio is big, as in very small 
particles where Abuik > A sur f., the over-broadening should occur on the LDOS distribution 
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Figure 3: Knight shift spectrum of bulk sites when U = IV (thick line), compared with the 
case with U = (thin line). M = 7 x 7 x 7, M surf . = 49 and M bu i k = 294. W = IV and 
rj = 0.05V. 



of bulk sites. The cluster configuration must be such that the surface sites are in higher 
proportion with respect to the bulk sites. This situation can be obtained for a structure with 
6x6x6 orbitals, taking all faces of the cube with shifted site energies U = 2V (M sur f. = 152 
and Mb u ik = 64). In this case A bu ik = 3.5A surf .. Figure 4 shows the Knight shift spectrum 
of bulk sites for this new configuration. The distribution (thick line) is broader than the 
corresponding homogeneous case (thin line). The reason is the same as the one discussed 
in connection with Fig. 1 provided that the relation A sur f./Ab u ik is now inverted. However, 
the NMR spectrum of surface sites (Fig. 5) exhibits a considerable over-broadening, unlike 
the case presented in Fig. 3. This is due to the fact that, because of the topology of surface 
sites, completely surrounding bulk sites, there is a strong mixing between surface states and 
interior states, which also produces two fluctuation scales on the LDOS distribution at the 
surface. 
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Figure 4: Knight shift spectrum of bulk sites setting U = 2V (thick line), compared with 
the case with U — (thin line). M = 6 x 6 x 6, M sur f. = 152 and (corresponding to the six 
sides of the cube) M bulk = 64. W = IV and tj = 0.05V. 



We now discuss the surface effects on the scaling law. Fig. 6 shows the relative Knight 
shift as a function of the scaling parameter a = 77/ A. The solid lines correspond to our 
results for particles of various sizes with surface energies on all faces of the cube. These 
curves correspond to particles of 7 x 7 x 7, 10 x 10 x 10 and 12 x 12 x 12 sites. Experimental 
results by Kobayashi jT3] et al on aluminum nanoparticles are also shown. The qualitative 
dependency of the relative knight shift with a agrees very well with our result. In addition, 
a deviation of the scaling law is clearly observed in both, simulations and experiment. Large 
particles lie further to the right than small particles. This can be interpreted as follow: due 
to the existence of surface states, the mean spacing level at the Fermi energy is no longer A, 
since this last spacing is defined according to the density of states of the bulk. Therefore, 
a is no longer the scaling parameter. If we define A' as the mean spacing level at the 
Fermi energy in presence of surface states and a 1 the corresponding scaling parameter, then, 
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Figure 5: Knight shift spectrum of surface sites when U = IV (thick line), compared with 
the case with U = (thin line). M = 6 x 6 x 6, M surf . = 152 and M bulk = 64. W = IV and 
r] = 0.05V. 

according to Eq. El a = a' 4^. Since bulk states and surface states become uncorrelated and 
consequently the repulsion of levels around the Fermi energy decreases, A' < A. Therefore, 
the factor ^ counts for the deviation of the scaling law. It is smaller for smaller particles 
than for larger particles which is consistent with the experimental and predicted results in 
Fig. 6. 

4 Conclusions 

According to our results, we expect that the NMR spectra of either surface or bulk nuclei in 
small metallic nanoparticles, manifest the existence of surface states as an over-broadened 
line which could also be interpreted as a strong disorder. In fact, the deviation of the scaling 
law shown in the experiment on aluminum nanoparticles ^H] is an indication of different 
distributions of energy spacings coming from bulk and surface states. Additionally to the 
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Figure 6: Relative Knight shift x max . vs. the scaling parameter a. The solid lines correspond 
to simulations of various particle sizes with surface energy U = 2V for all faces of the cube. 
Experimental values correspond to aluminum nanoparticles |15j . In both simulations and 
experiments a deviation of the universal scaling law is observed. 

over-broadening of crystalline Ni particles [Zj which is consistent with this view, there are 
other ways in which these effects might be observed and tested: the NMR technique SEDOR 
(spin echo double resonance) consists of simultaneous irradiation of resonances of nuclei of 
two different species which are coupled either by direct dipolar interaction or through the 
conduction electrons. With this method the contribution to the resonance spectra of the 
nuclei from the surface can be inferred. A very interesting application ^2] is the test of 
Platinum nanoparticles using Carbon as a local probe. The last is chemisorbed on the 
Pt surface as CO. There, the presence of a second peak in the Pt SEDOR line might be 
interpreted as the coexistence of surface and bulk states. 

Also based on the CO chemisorption at the surface of a metallic nanoparticle, it has been 
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observed [T71 ITS! IT9*] that Carbon NMR presents a shift toward a metallic nature. Not only 
the Korringa equation for the relaxation time T\ is verified (1/Ti oc (N(ep,r)) 2 /kBT) but 
the line shape presents the inhomogeneity proper of quantum size effects. That is, due to 
the mixing between the conduction band of the metal and the CO molecular orbitals, the 
frequency shift of the C nuclei contains information of the LDOS on the surface nuclei . 
Eventually an NMR study of chemisorbed molecules also would give evidence of anomalous 
fluctuations in surface LDOS. 
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